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1. Introduction. 

The in-depth understanding of laser-plasma interaction is of paramount importance 
for the eventual success of inertial confinement fusion research, but is also interesting 
for magnetic confinement fusion research, since tokamak plasmas can be heated by 
electromagnetic waves. The complex kinetic phenomena involved in this interaction, 
and the instabilities they may generate P need to be studied by kinetic models 
even though hydrodynamic models |3] are more affordable to simulate complex, high- 
dimensional geometries. However, the use of the full 3D Vlasov-Maxwell system is 
of course impossible in most practical situations. Therefore, the reduced Vlasov- 
Maxwell system for laser-plasma interaction (hereafter called the "laser-plasma system"; 
see (PHHl) below) was introduced in j2]. The model has been shown to capture some 
essential features of this interaction [21 C] ; and it has been successfully used for deriving 
relevant physical models in novel situations |3]. 
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The laser-plasma system has been the object of several mathematical investiga- 
tions El ■ In this framework, it is interesting to find classes of exact solutions which 
may serve as "reference solutions" , to which other solutions may be compared in order 
to study the dynamic of the interaction. Reference solutions for the Vlasov-Poisson 
system are, for instance, the Bernstein-Greene-Kruskal or BGK modes jSj, given by a 
distribution function of the form fiW), where W is the energy of one particle. When the 
function / is convex, as in the Maxwellian case f{W) oc e~^/^, such solutions represent 
fundamental equilibrium states; their existence and stability are well known [9 . For a 
general /, BGK solutions may represent various wave phenomena; they have been the 
object of many investigations in the physical community PHl and references therein]. 
The mathematical theory is still less developed; interesting existence and (in)stability 
results have appeared recently ^21 ^] • For both Vlasov-Poisson and Vlasov-Maxwell 
systems, there are also hnearised solutions leading to the dispersion relations of the 
various types of waves; e.g. for electromagnetic waves u"^ = uj'^ + fc^, where (cu, /c) are 
the pulsation and wave number, and cUp is the plasma pulsation. 

In this article, we shall introduce a class of exact solutions to the laser- 
plasma system which generalises, at the same time, Maxwellian equilibria and linear 
electromagnetic waves. Indeed, we investigate the existence of quasi-static solutions 
where the distribution function is at any time proportional to the Boltzmann factor; 
this static character can be reconciled with the electromagnetic character of the system 
by assuming a harmonic time dependence of the electromagnetic field and a circular 
polarisation. This ansatz was already used in but in a different physical and 
mathematical context. The latter work investigates the existence of solitons in an 
electron-positron plasma, where no charge separation occurs. Here we are dealing with 
a general ion-electron plasma, and we are looking for space periodic solutions. 

The paper is organized as follows. We recall the mathematical results known about 
the laser-plasma system and introduce the quasi-static model in section El Then, in 
sectionjniwe solve (in the space periodic setting) the so-called Boltzmann problem, which 
consists in finding the equilibrium density given the electromagnetic potentials, and we 
estimate its solutions. In section jH we construct a fixed point application and we study 
its properties. The existence of Boltzmaniann equilibria then follows by applying the 
Schauder fixed point theorem. Several extensions of the model are briefly discussed in 
section O and we conclude in section (HI 

2. The harmonic Boltzmannian model. 

The reduced Vlasov-Maxwell system for laser-plasma interaction describes the evolution 
of the distribution function of a population of electrons in a one space dimensional 
plasma interacting with a laser wave. In a first approach, we assume that the ions are 
at rest and their density is given — which is physically acceptable at the time scale of a 
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laser wave. After a suitable rescaling [3], this system can be cast in the following form: 
df P df / , , A{t,x) dA\ df 

dE , , , . dE .. , , . 

-Q^ = Ph{x) - p{t,x), — -j(t,x) = 0, (2) 

d'^A d'^A 

^-^ + P(^'^)^(^'^) = 0' (3) 

where: f(t,x,p) is the electron distribution function {p denotes the x-component of the 
momentum vector); E is the x-component if the electric field; A = {0,Ay,Az) is the 
vector potential of the laser wave; ph{x) is the (static) background ion density; 71,72 
are Lorentz factors. We distinguish three cases: 

(i) the non-relativistic case (NR), 71 = 72 = 1; 

(ii) the quasi-relativistic case (QR), 71 = (1 +p^)^''^,72 = 1; 

(iii) the fully relativistic case (FR), 7i = 72 = (1 +p2 + |A|2)i/2^ which is the original 
model of [2]. 

The moments p, p, j are given by 

p{t,x):= / f{t,x,p)dp, p{t,x):= / — — dp, j{t,x) := / —f{t,x,p)dp. (4) 

Jr Jr 72 Jr 7i 

We supplement the system dU |21 El) with initial conditions 

f{0,x,p) = fo{x,p), {x,p) iE,A,dtA){0,x) = iEo,Ao,A,), xeR. (5) 

In it was proved that, for suitable initial conditions, ((IHSl) has a unique classical 
solution, which is global in time in the QR case, and local in time in the NR case. In 
the latter case, the classical solution can be extended to a global weak solution with 
/ continuous and A continuously differentiable in all their variables. The FR model 
was studied in [6 . It was shown that admits a unique global classical solution 

preserving the total energy. The stationary solutions of these models in a bounded 
domain have been analysed in [7]. 

All three models admit space periodic solutions. If the initial data are L-periodic 
in X and satisfy the neutrality condition 

/ / fo{x,p) dpdx = / p^{x) dx =: M, (6) 
Jo Jr Jo 

then, by using the continuity equation dtp + dxj = 0, we deduce that the system remains 

globally neutral at any time t > 

pL p pL p pL 

/ / f{t,x,p) dpdx = / / fQ{x,p) dp dx = / py,{x) dx. (7) 
Jo Jr Jo Jr Jo 

By uniqueness of the solution one gets also that {f{t),E{t),A{t)) are L-periodic in 

space for any t > 0. From now on, we work in the framework of periodic functions: 

all differential equations will be implicitly supplemented with L-periodic boundary 

conditions. 



Boltzmannian waves in laser-plasma interaction 



4 



From ((7j) we deduce the existence of a unique function V = V{t,x), satisfying 
dlV{t,x) = p],{x) -p{t,x), V{t,0) = and {V,d^V){t,x) = {V,d^V){t,x + L), for all 
{t,x) G [0, +oo) X R. The field E derives from the potential V, i.e., E = d^V. 

The purpose of this article is to study the existence of particular solutions of (HI |21 
IHl) corresponding to local Boltzmannian equilibria. These are defined by f{t,x,p) oc 
Q-w(t,x,p)/e -yyj^gj-g W{t,x,p) is the energy of one particle being at the phase space 
point {x,p) at time t, and 9 is the scaled temperature. As it is well known, such 
functions are solutions to the Vlasov equation ((T)) iff is independent of time. Thus, 
we assume that V does not depend on t, and that A is time-harmonic and circularly 
polarized, i.e., 

Ay{t, x) + i Azit, x) = a{x) e"^*, with a priori a{x) G C. 

Then the energy W{x,p) is given, according to the relativistic character, by 
W{x,p) = i (p2 + \a{x)\^) + V{x), in the NR case, 
W{x,p) = a/1 +p2 + i |a(x)p + V{x), in the QR case, 
W{x,p) = a/1 + + + V{x), in the FR case. 

Imposing the constraint ® yields 

Q-w{x,p)/e 

f{x,p) = M ,,,, , ^ V(x,p)GMl (8) 

By direct computation we check that in all three cases / solves the Vlasov equation (^. 
We then observe that j{x) = ^f{x,p) dp = 0, for x G M, and thus the system (jH 121 
El) reduces to 

V"{x) = pb{x) - p{x), xeR, (9) 

— co'^a(x) — a"{x) = —p{x)a{x), x G M. (10) 

with p = J^f dp, p = J^^ dp and / given by ©. 

Of course, we are interested in solutions such that a ^ 0, otherwise we find a 
Vlasov-Poisson equilibrium. If such a solution exists, a appears as an eigenfunction of 
the operator Ap := — ^ + p{x), associated to the eigenvalue tu^. It is well known that 
these eigenvalues are real and generically simple; in particular, the lowest eigenvalue 
is always simple. As the coefficients of Ap are real, we infer that both 3fJ(a) and 
9(a) are eigenf unctions; thus, generically, they must be proportional. In other words, 
a{x) = a(a;)e"^, where a is a real eigenfunction and ^9 G M. Then, |a(a;)p = a(a;)^, and 
W, f, p, p are defined in terms of a; while we may take (p = hj rotating the axes 
Oy, Oz. This means that, without loss of generality, we may restrict our search to real 
functions a solution to (|TO|l . 

We now rewrite the model Q ITn|) in a form which will prove more convenient for 
analysis. We shall denote by the subscript # the spaces of L-periodic functions, e.g.: 
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■■= {g e Ll^{R) : Vx, g{x + L) = g{x)}, (R) := {w G ^"(R) : Vx, w{x + L) = 
w{x)}. First, we introduce the operator $ : i^^(R) —>■ C^(R) given by 

<^[g]=w eCl{R), -w"{x) = g{x), xe{0,L), w{0) = w{L) = 0, 

for any g G L^(R). Then, we consider the function : R ^ R given by 

e-^(^)/^ = / exp _ W(-^P\-n^) dp, X G R, 

namely, according to the relativistic character: 

NR: e-^^^^/' = ^-^^-f I''' [ e^^'/^^ dp, (11) 



QR: e-'^(^)/'^ = e-'^(^)^/2'^ / exp-^^^i^ dp, (12) 



FR: e-^(^)/^= / exp - ^il^l^ dp. (13) 

Notice that there is a constant C{6) G R such that 

^(x) = ^^1^ + (:7(^), in the NR and QR cases. (14) 
In the FR case, by observing that 



2 

we obtain 



i(v/l+p2 + |a(x)|) < Vl+p2 + a(a;)2 < v^rTp2 + |a(x)|, 
C2(^) e-l'^(^)l/^ < e-'^(^)/^ < Ci(^) e-l''^^)!/^^ 



with 



Ci(0) := / exp-^^^^2dp> / exp-^/^dp=:C2(^). 
Jm -^(7 



Finally one gets 

^ < i;{x) + e\nC^i9) < \a{x)\ (15) 

The density p can be expressed in function of ip; and the system (jS HTUI) can be 
recast as: 

/(x,p) = ire-^("'^')/^ 

p{x) = Ke -0 , xgR, V^=$[p-pb], (16) 

p{x) = K [ le-^("'^')/' dp, (17) 
Jr 72 

—u'^a(x) — a"{x) = —p{x) a(x), (18) 

where the constant K = M (^J^ Q-w{y,q)/9 dg dy^ is such that p{x) dx = M. 
We call this system the Boltzmann-Helmholtz equations; they can be seen as a sort of 
non-linear eigenvalue problem. 
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3. The Boltzmann problem. 

For the moment we suppose that the function if) is given and we solve the so-called 
Boltzmann problem (fTBj) . The proof of the following proposition is immediate and left 
to the reader. 

Proposition 1 For any function g G L^(]R) we have: 

||<^'MI|l-(r) < L\\g\\Li{o,L)- 
If the function g satisfies g{x) dx = 0, then ^[g] G C^(IR) and we have: 



< 



\9\\l^o,l)- 



Proposition 2 Let ip e L#(M), G Mb > 0, M = Uk{x) dx and 6 > 0. 

Then there is a unique function u G -^^^(IR) such that: 

(19) 



u = M 



Moreover it satisfies: 



exp 



!o exp 



■>p(y)+^u~u^]{y) 



dy 



0<M<mf— exp-U {i,{y)-C)dy-M{^-C)+ALM\=:u^- (20) 



and zf^pe iyi'°°(M) then u G iy^'°°(M) and: 

Upip + 2M 

Lip u < u^. 

U 



(21) 



Proof: One readily checks that (jl9p is equivalent to the minimization of the functional 



J[v] 



{9a{v{x)) + 



dx 



<^[V- Ub] 



+ ip{x)v{x)} dx, 



under the constraint v{x) dx = M, where a{s) = sins, s > and o"(0) = 0. This 
problem is a variant of that considered in ^] and its well-posedness follows from a 
similar argument. The functional J is strictly convex, l.s.c. and bounded from below 
on the set ^ 

]C{L, M) = jt; G L^(R) : v>0, v{x) dx = m| . 

Indeed, by applying the Jensen inequality with the convex function cr, the measure 

d;x = e-^(^)l() ( e-^^y)ie dx and the function v/e'^'/^, one gets: 



J[v]> / {ea{v{x)) +i){x)v{x)} dx>eM\n 



M 



dy 







saying that inf„g^(j;,^M) J[v\ > — oo. Take a minimising sequence («„)„. By using the 
Dunford-Pettis criterion we can assume (after a suitable extraction) that {un)n converges 
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weakly in L^(0, L) towards a function u G 1C{L, M). Since J is convex we can pass to the 
limit by involving the semi-continuity of J and we obtain that J[u\ = in{^^ic(L,M) <^[^]- 
Writing the Euler-Lagrange equation we obtain 

6(1 + Inw) + - Mb] + ^ - a = 0, 

where a enters as the Lagrange multiplier associated to the constraint u{x) dx = M 
and thus we deduce By using now the Jensen inequality with the convex function 
1 1— >■ e~*, the measure dfi = dx and the function {ip + $[m — Ub]) /6 we obtain: 



exp 



L 



+ - Ub] 

e 



dy 



< j exp 
'o 



ip + ^[u — U]^]\ dy 

T' 
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Therefore by using Proposition ^ we deduce 



e dy 



1 1 



_i£+«K-_uiJ 1 11 

e « d?/<-exp-(- 



ij dy + 2LM), 



which implies 



M 1 /I /"^ \ 

< u{x) < — exp - y ilj{y) dy - ^(x) + 4 L Mj . 

This inequality is unchanged by replacing with — C, for any C G M; one thus 
infers fl20|) . Assume now that ^/^ G iy^'°°(]R). By taking the derivative with respect to x 
in (imjl one gets by using Proposition [T] 

and (PT|) follows immediately. □ 



4. The fixed point application. 

For any a > we define the fixed point application JF^ : iy^'°°(]R) 
for any a G W^'^iR) where: 

• T/^ is given, according to the case, by (fTTj) . (fT^ or (fT^: 

• p is the unique solution to the Boltzmann problem 

rL 

Vi+$[p-Pfa] / 

p = Ke 9 , / p{x) dx = M ; 
Jo 



p = p in the NR and QR cases, while in the FR case: 
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• A is the first eigenvalue of tlie operator Ap = — ^ + p with L periodic boundary 
conditions, i.e., 



X = inf JO " y ' ■ ry ^ . ^22) 



• a is the corresponding eigenfunction of Ap-. 

-a"{x)+p{x)a{x) = Xa{x), xG(0,L), (23) 

a(o) = a(L), a'(o) = (24) 

such that a > and a(a;)^ da; = a^. 

Remark 1 It is well known that the first eigenvalue of Ap with L-periodic boundary 
conditions is simple USf and that the eigenfunction vanishes nowhere. Therefore 
a = J-'att is well defined. 

The properties of the apphcation JF„ are summarized up below. 

Proposition 3 Assume that pb ^ -^#(1^); Pb > 0, Jq Ph{x) dx = M and let a, 6 be 
positive real numbers. For any a G W^°°{Wj such that a{xy dx < construct 
ip, p, p, A and a = J-'aa as above. 

(i) There are constants Pi,,a^ depending on a,L,M,6 such that 

ii~ii / II II / II /|i ^ ||a||L°°(R)||a'IU=°(R) + 2M 
||P||l°°(r) S ||P||l°°(r) S P*, ||P IU°=(r) S ^ P*, 

ll~/|| , ||a||L°°(M)||o.'||L°°(R)(l + ^) + 2M n^\^ ll~ll / 

IIP ||l°°(r) S 7 Pi., U S A S P*, ||a||w2,^(]R) < a^. 

u 

(a) J-'a is continuous with respect to the topology o/C^(]R) on the set C = {a G C^(IR) : 
lkllL2(o,L) < a and ||a||i,oo(R) + ||a'||2,cx)(iR) < a^}. 

Proof: (i) Take a G W^°°{M.) such that ||a||i2(o,L) < a. In the NR and QR cases, we 
deduce from (fT^ and ^lU^ the bound 

< P < ^ exp- (- I ^ dx + 4LMj < - exp- (-a^ + 4LMj . 

In the FR case, combining ()15|) and ()20|) yields 

< p < — exp - f - y {i){y) + e\nCi{d)) dy - m.i{i) + d\YiCi{e)) + AL M \ 

We check easily that in all three cases we have l^/^'l^^)! ^ a; G M and thus, 

by Proposition 121 we deduce 

II /|| ^ ||a||L°°(R)||a'||Loo(]j) + 2 M 

IIP ||l°°{r) ^ J, ||P||l°°(r)- 
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The estimate for p follows, since in all three cases 72 > 1 and < p{x) < p{x). By taking 
the derivative with respect to x in the expression of p we obtain by direct computation 

,. / l\ 2M\ 

\\P ||l°=(r) < I ||a.||L°°(K) ||a ||l°°(ir) \ ~^ ) ~^ 0~ ) 

We now estimate the eigenvalue A and the eigenfunction a. Equation (j221) shows that 
A > and, by taking 6 = 1, A < p*. Then, from we deduce 



L pL 

{a'{xY + p{x) a{xY} dx = X a(x)^ dx, 
Jo 



and hence ||a||_ffi(o,L) < VA+T a. By using the Sobolev inclusion H^{0,L) C L°°(0,L) 
one gets easily that ||a||Loo(R) + ||a'||Loo(K) + ||a"||Loo(K) < a^{a, L, M,9). 

(ii) Take a sequence (a"')„ C C which converges towards a E C with respect to 
the topology of C^(]R). For any n let ip"", p", p". A", a" = jF^a" constructed as in 
the definition of the fixed point application. Similarly let ip, P) p. A, a = Taa- The 
sequence {ip^)n is bounded in iy^'°°(M) and therefore, by the Arzela-Ascoli theorem 
we can extract a subsequence converging in C^(M). Obviously the limit function is ip 
and by the uniqueness of the limit we deduce that the whole sequence (V'")n converges 
towards ip in C^(M). In the same manner, since sup^dlp^llicx))-]^) + || ^p"||Loo(]g)) < +00 
we deduce that p^ ^ p, p^ ^ p in C|^(M). 

The fact that lim„^+oo A" = A stems from general spectrum continuity 
theorems [Hj, or can be directly deduced from (j22|l . Finally, as sup„ Ha^H vi/2.°°(r) < +00, 
we can extract a subsequence {a^'')k converging in C^(M) towards some function a. By 
passing to the limit with respect to k in the weak formulation of a"* we obtain that the 
limit a satisfies 

—a (x) + p{x) d{x) = A a(x), x G (0, L), a(0) = a{L), a (0) = a (L). 

Moreover since a" > 0, ||a"||2,2(o,L) = ci for any n, we have a > and ||a||L2(o,L) = a and 
thus a = a = J^aP'- By the uniqueness of the limit we have lim„^+ooa" = a in C^(R). 
□ 

We are now in position to prove our main result by using the fixed point method. 

Theorem 4 Assume that pb G L^(M), pb > 0, ph{x) dx = M and let 6 be a positive 
real number. For any a > there is at least one classical solution (p, a) G C^{R)xC^(M.) 
for the Boltzmann-Helmholtz equations ^^^^ satisfying p > 0, p{x) dx = M, 
a > 0, a{xY dx = a^. 

Proof: Consider JF^ = J-'a\c- The set C is convex and compact in C^(]R); by 
Proposition El we know that fiaiC) C C and that JF^ is continuous with respect to 
the topology of C^(M). By the Schauder fixed point theorem we deduce that there is 
a fixed point a G C. By construction we have a > 0, a{xy dx = a^. Consider now 
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ip, p, p, A as in the definition of jF^a. Obviously A > 0, p > 0, /q p(x) dx = M and 
we check easily that (p, a) G C^(]R) x C|t(M). Observe that A > 0. Indeed we have 

/p^{a'(x)^ + p(x) a(x)^} dx ^ p{x) a(x)^ da; 
a(x)^ dx a{xy dx 

If A = then p(x) a(x)^ = for any x, and since by construction p > we deduce that 
a = which contradicts a(x)^ dx = > 0. Consider now uj = \/X > and thus 
(p, a) is a solution of (fTBHTS|) . □ 

5. Extensions. 

One could investigate the existence of "non-linear harmonics" of the "fundamental 
mode" given by Theorem IH i.e., solutions to |T8|) where u'^ is not the first eigenvalue 
of Ap, but one of higher rank. Unfortunately, it appears impossible to generalise the 
construction of JF^ to these eigenvalues. The reason is that, with periodic boundary 
conditions (unlike the Dirichlet, Neumann or Fourier b.c), these eigenvalues may be 
double for some "exceptional" densities p. For instance, if pb = est and a = est, then 
p = est and all eigenvalues except the first one are double. There is apparently no 
way of defining a continuous mapping p i— > a in the neighbourhood of the exceptional 
densities. Nevertheless, the existence of harmonics is very likely, as the eigenvalues are 
generically simple. 

Another interesting extension is the case where the ion density is no longer given, 
but is also proportional to the Boltzmann factor. Let us denote by the subscript 1, 
resp. 2, the quantities relative to the electrons, resp. ions; we introduce a new 
parameter p representing the electron/ion mass ratio. Then, we have /i oc e~^^/^^ 
and /a oc e-^^^^\ The energy Wi of one electron is given in section |21 that of one ion 
is, according to the relativistic character: 

W2{x,p) = lpip^+ |a(x)n - Vix), (NR), 

W2{x,p) = p-' y/l + {ppy + lp\a{x)\^ - V{x), (QR), 

W2{x,p)=p-'^/lTp^iJ^TH^-V{x), (FR). 

We arrive at the following system: 

V"{x)=p2{x)-pi{x), xeR, (25) 

-uj'^a{x) - a"{x) = -{pi{x) + pp2{x))a{x), x eR. (26) 

The arguments of sections El and El can be extended without bad surprises to this two- 
species model. However, the solutions corresponding to the first eigenvalue are not 
very interesting: one easily checks that all the functions pi, p2, pi, P2, a are constant, 
and V = 0. 
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6. Concluding remarks. 

In this article, we have shown the existence of quasi-equilibrium solutions to the laser- 
plasma system, where the distribution function is Boltzmannian and the electromagnetic 
variables are time-harmonic, at least at the fundamental frequency. The existence of 
solutions at higher frequencies is probable, both for one-species and two-species models. 
These solutions appear as generalisations of Vlasov-Poisson equilibria, but are clearly 
different from them as an electromagnetic wave is present. They can be viewed as a 
simple case of non-linear interaction between the electron plasma oscillations and the 
laser wave. The implicit relation (through the spectrum of the operator Ap) between 
the frequency uj and the space period L yields in the linear limit the dispersion relation 
for electromagnetic waves. 

Quasi-equilibria can serve as references for analysing the dynamics of laser-plasma 
interaction, e.g. Raman and Brillouin scattering, which are among the most challenging 
issues to deal with in order to achieve controlled inertial confinement fusion. Indeed, 
from a dynamical point of view, it should be noted that quasi-equilibria may be unstable, 
unlike the Vlasov-Poisson equilibria which are non-linearly stable, even under ID 
Vlasov-Maxwell perturbations [5 . These solutions may also serve as benchmarks for 
testing numerical codes, even though the numerical solution of the Boltzmann problem 
appears quite difficult when Mb and/or ifj feature large variations. 
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